The Immunity of Polymer-Microemulsion Networks 



Guy Hed and S. A. Safran 
Department of Materials and Interfaces, Weizmann Institute of Science, Rehovot 76100, Israel 

The concept of network immunity, i.e., the robustness of the network connectivity after a random 
deletion of edges or vertices, has been investigated in biological or communication networks. We 
apply this concept to a self-assembling, physical network of microemulsion droplets connected by 
telechelic polymers, where more than one polymer can connect a pair of droplets. The gel phase 
of this system has higher immunity if it is more likely to survive (i.e., maintain a macroscopic, 
connected component) when some of the polymers are randomly degraded. We consider the dis- 
tribution p(o) of the number of polymers between a pair of droplets, and show that gel immunity 
decreases as the variance of p(o~) increases. Repulsive interactions between the polymers decrease 
the variance, while attractive interactions increase the variance, and may result in a bimodal p(cr). 



I. INTRODUCTION 

The immunity of networks refers to the ability of a 
network to preserve its function, even when a certain 
number of edges or vertices is removed. This attribute is 
of interest in many types of networks . The immunity 
of genetic, cellular networks reflects the robustness of the 
network under a mutation or deletion of one of the genes 
0. The robustness of the metabolic network of single 
cell organisms is demonstrated by their persistence and 
growth despite environmental interventions 0. In com- 
munication networks, immunity refers to the ability of 
the network to pass messages among most of its vertices 
even after some communication lines are cut, or some ver- 
tices are damaged 0. A communication network is im- 
mune to such an attack, if it retains a connected section 
that spans enough vertices to allow some communication 
(even if quite indirect) between any two vertices. 

In this paper, we extend these ideas to analyze the im- 
munity of a network composed of droplets connected by 
polymers . The gel phase of such a system is character- 
ized by a connected part of the network that spans the en- 
tire volume. The immunity of the gel is characterized by 
the fraction of polymers that can be removed (by chem- 
ical breaking, depolymcrization or by shear forces) from 
the network while the gel maintains its macroscopic char- 
acteristics, and does not transform into the fluid phase. 
That is, a gel is immune if the physical network that 
comprises the gel retains a macroscopically connected 
component even though a relatively high fraction of the 
polymers that form the network have been deleted. 

A simple polymer gel is a network of connected poly- 
mers in which each edge (i.e., bond) of the network con- 
sists of one polymer; in this case, the concept of immu- 
nity is idential with that of percolation 0. However, in 
the system of polymers that join microemulsion droplets, 
each edge (i.e., bond) that connects two oil droplets may 
contain more than one polymer - and sometimes many 
polymers. The percolation problem deals with the re- 
moval of entire edges from the network, while immunity 
refers to the removal of individual polymers from a given 
edge in systems where there can be many polymers in an 
edge. 



Examples of equilibrium, network-forming systems in- 
clude surfactant solutions, gels of biological molecules or 
synthetic polymers. A particularly elegant experimental 
realization of a transient, self- assembling, physical net- 
work has been reported in . The system consists of oil- 
in-water microemulsion droplets connected by telechelic 
polymers; the latter have a hydrophilic backbone with a 
hydrophobic group at each chain end. The insolubility 
of these groups in the continuous water phase and their 
solubility in the oil droplets results, in some cases, in a 
network of droplets connected by the polymers. Mixtures 
of telechelic polymers and microemulsions have a wide 
range of technological applications, including paints, cos- 
metics and enhanced oil recovery. Precise control of the 
structural and rheological properties of such materials is 
essential for effective and reliable performance. Such con- 
trol is possible using the telechelic additives, that form a 
transient network with controlled rheological and struc- 
tural properties that depends only on the physical prop- 
erties of the system such as the size and concentration of 
the polymers and the droplets . 

Apart from its applied interest, the telechelic- 
microemulsion mixtures serve as model systems for the 
understanding of a more general class of transient net- 
works. The advantage of this particular system is that 
the parameters that control the thermodynamics and 
structure can be easily identified and independently 
controlled: the concentration of possible vertices (the 
droplets) and the connectivity of the network (the num- 
ber of polymers per droplet). 

The experimental systems exhibit a phase transi- 
tion from a single phase to coexisting high density and 
low density network phases. Zilman et al. showed 
that this transition is expected from entropic considera- 
tions alone, and occurs even without any particular ener- 
getic interactions between the droplets and polymers. In 
this paper, we extend the analysis of and examine the 
immunity of such polymer gels to random degradation of 
the polymers such as those that can result from chemi- 
cal degradation or from shear forces. As the polymers are 
removed, the macroscopic network becomes disconnected 
and the gel reverts to a liquid. We show the interactions 
between the hydrophobic groups at the ends of the poly- 
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mers have an important affect on the immunity and the 
phase behavior of the system. 

The relation between the degradation of the polymers 
and the breakdown of the network is not trivial. Each 
edge in the network may contain several polymers, each 
of which connect one pair of droplets. The fraction of 
polymers that must be degraded in order for the network 
to lose its connectivity depends on the distribution, p(er), 
of the number of polymers in an edge. As the variance 
of the distribution increases, the network becomes more 
sensitive to polymer degradation, and will lose its connec- 
tivity when a smaller fraction of polymers is degraded. 

We consider the distribution, p(a), for various physi- 
cal situations. Since we are treating relatively rigid and 
short polymers, we do not consider the interactions be- 
tween the chains, but focus on the interactions between 
hydrophobic groups at the end of the polymers, which 
we denote as stickers. In particular, we predict the ef- 
fect of both repulsive and attractive interactions between 
the stickers within a droplet. Repulsive interactions are 
expected from entropic and excluded volume consider- 
ations, while attractive interactions can arise from Van 
der Waals forces, or from the local change of curvature 
that the hydrophobic ends induce on the droplet surface. 
We find that repulsive interactions decrease the variance 
of p{u) while attractive interaction increase the variance. 

Attractive interactions between the stickers may lead 
to the co-existence of dense edges, with high density of 
polymers, with dilute edges, with low polymer density. In 
that case, the distribution of the edge occupation num- 
ber, p(a), is bimodal, and its width is of the order of its 
average. As the system approaches the critical point as- 
sociated with the emergence of this two-phase behavior, 
the fluctuations in the number of polymers in an edge 
as well as the variance of p(cr), increase. Thus, as we in- 
crease the strength of the interaction between the stickers 
toward its critical value, local concentration fluctuations 
increase, and the immunity decreases. 

In Sec. |H]wc describe the underlying physical model 
for the network forming systems of interest, and the sim- 
plified theoretical model we use in our calculations. Next, 
in Sec. IIIII we examine some simple distributions ior p(a) 
and demonstrate the dependence of the immunity of the 
gel on the variance of the distribution. The effects of re- 
pulsive and attractive interactions between the polymer 
ends on or in a droplet are presented in Sec. IIVI We use 
an expression for the Landau free energy of the system 
of both interacting polymers and droplets to obtain the 
effect of the interactions on the variance of p(cr), and we 
predict the phase diagram of this system. 



II. THE MODEL 

We map the physical space to a discrete regular lattice 
, where each site can be occupied by an oil droplet. 
Each pair of neighboring droplets may be connected by 
polymers. The set of polymers that connect such a pair 



is referred to as an edge of the network, whose vertices 
are the droplets. 

From the applications of percolation theory to gelation 
we know that if the fraction, k, of occupied edges exceeds 
the percolation threshold, k c , of the lattice, the system 
is in the gel phase. In the usual discussions of gel con- 
nectivity and percolation, each edge is occupied by one 
polymer. After the degradation of a fraction q of the 
polymers, the condition for the survival of the gel phase 
is k(1 — q) > k c - 

However, the system we discuss here has two impor- 
tant and interesting differences from the usual percola- 
tion models: 

1. The vertices of the polymer network only exist at 
those lattice sites at which the oil droplets are 
present, so only two neighboring vertices of the lat- 
tice that are each occupied by droplets can be con- 
nected by an edge. In the following, we refer to such 
a pair as an edge even it is not connected by poly- 
mers, in which case we consider it as a disconnected 
edge. If a given droplet has a neighboring vacant 
lattice site, we refer to the the droplet surface that 
face this site as a free face. 

2. An edge may be occupied by a number, a, of poly- 
mers, whose value is determined by volume lim- 
itations (in the case of excluded volume interac- 
tions), explicit interactions between the polymers 
or between the stickers, and the polymer chemi- 
cal potential - which is fixed by their concentration 
in the solution. The polymer chemical potential 
determines the average number of polymers in an 
edge, (cr). 



III. NETWORK IMMUNITY AND EDGE 
OCCUPATION DISTRIBUTION 

In this section, we concentrate on the effect of the 
edge occupation distribution, p(<r), on the immunity. We 
do not consider here the droplet-related degrees of free- 
dom, but rather assume that they are fixed in a given 
(static) configuration. The following section considers 
the droplet-related degrees of freedom. The only degrees 
of freedom we relate to in this section are occupation, a, 
of the edges between neighboring droplets. 

If the droplets are dense enough, and most of the edges 
arc connected by one or more polymers, then the system 
has a connected component that spans its entire volume, 
and the system is in the gel phase. If the fraction of dis- 
connected edges (i.e. those that contain no polymer), is 
large, the system is in the fluid phase. A system in the 
gel phase that undergoes a degradation of a fraction of its 
polymers, has higher immunity if fewer edges are com- 
pletely disconnected by the degradation process. This 
will depend on the number of polymers per edge and on 
the variation of this number throughout the network. 
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For example, if the system has a distribution of edges 
in which some edges have large numbers of polymers and 
some have very few, it will have lower immunity than a 
system in which the number of polymers in each edge is 
the same. This is because the random removal of poly- 
mers from edges with lower than average occupancy more 
readily leads to the complete disconnection of that edge, 
as demonstrated in Fig. ^ 

The probability for a given edge to be occupied by 
<7 polymers is p(cr); the fraction of unoccupied, discon- 
nected, edges is given by p(0) - the probability of a pair of 
neighboring droplets to have no polymers linking them. 

After a fraction q of the polymers is degraded, the new 
occupation of an edge which had a polymers is propor- 
tional to the product of the probability of a number, a' 
of polymers to remain in the edge, (1 — q) a and the prob- 
ability of the remaining polymers to be removed, q( a ~ a \ 
The resulting distribution can be written 

a— a' ^ ^ 

The new fraction of unoccupied edges is p(0) = 
'^2 (T p(&)<l' 7 ■ The higher this quantity for a given initial 
average edge occupation, (er), and a given degradation 
probability, q, the lower the immunity of the system. 

A. Simple distributions 




FIG. 1: Models of polymer-microemulsion systems with differ- 
ent distributions, p(cr), of the number of polymers in an edge. 
We compare a homogeneous system (a) with 11 polymers in 
each edge (mean m = 11 and distribution width s — 0) vs. a 
system with bimodal p((j) (b), as given by Eq. [I] (rn — 11, 
s = 9). After a random deletion of a fraction, q = 0.85, of 
the polymers, most of the edges in the homogeneous system 
(c) survive, and it still retains a macroscopically connected 
component. In the inhomogeneous system, most of the dilute 
edges are eliminated by the deletion (d) and the system no 
longer has a macroscopically connected component. 



In order to elucidate relation between the variance of 
the edge occupation distribution p(a) and the immunity 
of the gel, we first study two simple distributions. The 
first is a uniform distribution, 



Pip) 



2s + l 





m — s < a < m 
otherwise 



(2) 



with an average m and variance S% = (s + s)/3. 

If we randomly delete a fraction q of the polymers, then 
the fraction of unoccupied edges is 



p(0) 



1 



,2s+l 



2s + 1 



1 



(3) 



which is an increasing function of s for any value of q. 
Thus, in a distribution with a larger variance, the prob- 
ability of an edge to remain connected is smaller, and it 
is less immune to random polymer degradation. 

The second distribution we would like to consider is 
the bimodal distribution depicted in Fig. ^ 



p(a) = 



1/2 




fj = m i s 
otherwise 



(4) 



with an average m and variance s 2 . After a deletion of 
a fraction q of the polymers, the fraction of disconnected 
edges is 



p(0) = q m cosh(slog<7) 



(5) 



As the variance s increases, the system become less im- 
mune to random polymer degradation, as demonstrated 
in Fig. ^ This is relevant to polymer-microemulsion sys- 
tems with a bimodal distribution of the number of poly- 
mers in an edge, as expected when the attractive inter- 
actions between the hydrophobic stickers at the polymer 
ends are large (see Sec. II V|) . 

In the following sections we consider several physical 
models corresponding to different physical interactions 
between the polymers and stickers, that yield different 
edge occupation distributions, p{cr). We calculate the 
variance for each situation to estimate the gel immunity 
and how it depends on the interactions. 



B. Non-interacting polymers 

We first examine a naive, but instructive, model where 
any edge can accommodate an infinite number of poly- 
mers. In order to compare it with other models, we ex- 
press the variance as a function of the polymer concen- 
tration, which is the control parameter in most of the 
experiments First, we calculate the distribution in 
the grand canonical ensemble, in which the chemical po- 
tential, fi p , of the polymers is fixed. We then calculate 
the dependence of fi p on the polymer concentration, to 
determine fi p . 
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The distribution is given by 
p(a) = e~ e " P 



-<T/Jp 



(6) 



where the factorial in the denominator comes from the 
fact that the polymers are indistinguishable. This distri- 
bution has an average and a variance (a) = S% = e _Mp . 



C. Excluded volume interactions 

In a more realistic model, the excluded volume inter- 
actions limit the maximum number of polymers that can 
reside in a given droplet. There is thus a finite, maxi- 
mal number of polymers, ctm, associated with one edge. 
The volume restriction may arise from Helfrich-like, ex- 
cluded volume interactions between the polymer chains, 
or from the finite area occupied by each sticker in or on 
the surface of a droplet. 

If the only interactions between the polymers and be- 
tween the stickers are short range excluded volume in- 
teractions, we can assume that there are <jm sites in an 
edge, that can each be occupied by one polymer. The 
distribution of the edge occupation is 



p(a) = Z~ 



ctm\ - 



-cr/ip 



(7) 



where Z = (1 + e~ Mp )' TJU is the partition function. 

The free energy per site (we take all energies in units of 
ksT) is / c = — \og(Z)/<jM. The average edge occupation 
is 



{a} = a M 



<9/i F 



1 + e^p 



(8) 



and the variance of the distribution p(a) is 

d 2 fc <tm f, (a) 

cm 



<9/i p 2 



4 cosh (/ip/2) 



(9) 

The variance here is smaller than in the case of non- 
interacting polymers (Eq. 0). In the non- interacting 
case, the variance was equal to (a), where in the present 
case it is smaller by a factor of 1 — (a) joM- As expected, 
the difference between the distributions vanishes for the 
case om 3> (&)■ The narrowing of the distribution by 
the excluded volume interactions increases as the aver- 
age number of polymers per edge approaches aM, where 
the variance vanishes and the distribution approaches a 
delta function as a — > aM- This suggests that repulsive 
interactions between the polymers or the stickers increase 
the immunity of the gel under a random degradation of 
the polymers. 



IV. INTERACTING STICKERS 

We now calculate the free energy of a system of 
droplets and polymers in which the polymers show at- 
tractive or repulsive interactions in the multi-canonical 



ensemble, where the number of polymers and droplets is 
not fixed, and is determined by the chemical potentials 
/i p and /id, respectively. We shall consider the possibil- 
ity of polymers that connect neighboring droplets as well 
as polymers that form loops, in which both sticker ends 
reside in the same droplet. 

The physical model can be associated with the vertices 
(droplets) and edges of an underlying graph. The vertices 
are given by the sites of the underlying lattice which are 
occupied by the droplets. The edges of the graph con- 
nect pairs of neighboring occupied vertices on the lattice. 
Since each site may be occupied by a droplet, the statisti- 
cal ensemble we consider here includes the sub-graphs of 
the lattice (though most of them have negligible weight). 

We divide each droplet into faces, one for each edge 
of the graph that emerges from its vertex. If another 
droplet occupies a vertex that is connected to the same 
opposite edge, then the corresponding face is associated 
with the edge between the two droplets, and we call it 
a linked face. If there is no neighboring droplet at the 
vertex connected to a given edge, we term it a free face. 

We consider systems in which the radius of gyration of 
a polymer is smaller than the radius of a droplet. Thus, 
we treat each edge and each free face as a separate sys- 
tem. In an experiment, where the total number of poly- 
mers is fixed, the edges and faces interact through this 
constraint. In the multi-canonical ensemble, however, 
they do not interact. We can thus treat each edge or free 
face as an independent system and write its partition 
function explicitly. 

As discussed above, excluded volume interactions im- 
ply that a face can contain a maximum of aM stickers. 
This means that an edge has a maximal occupation of 
ctm polymers. These polymers can be links - with the 
two telechclic stickers of a given polymer attached to the 
two opposing faces of adjacent droplets, or loops - with 
both stickers of a given polymers attached to the same 
face of a single droplet. A face can contain a maximum 
of ctm/2 such loops. 



A. Polymer free energy 

We consider a model in which, except for the excluded 
volume interactions that determine the maximal num- 
ber of polymers per droplet, omi polymers interact only 
through their stickers (telechelic ends). The interactions 
between the polymer chains have a negligible contribu- 
tion for the case of relatively short and stiff chains that 
we consider here. Thus the energy difference between 
an edge that has two loops on the opposite faces of two 
neighboring droplets, and an edge in which these faces are 
linked by two polymers, is simply 2ei oop — 2eu n k, where 
ei 00 p is the configurational free energy of a loop, and eii n k 
is that of a link. We thus define an effective thermody- 
namic potential of a polymer in an edge 



M = M P - log(e~ 



(10) 
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An edge where all the polymers have a thermodynamic 
potential [i is equivalent to a one in which loops have a 
potential of \x v + e\ oop and links a thermodynamic poten- 
tial of /i p + eii n k, in the sense that both systems have the 
same partition function and hence yield the same distri- 
bution p(a). Thus, we assume that all polymers in the 
edge have a thermodynamic potential of /i, whether they 
are loops or links. 

In a free face, where the polymers are all in loop con- 
figurations, the thermodynamic potential of a polymer is 
M P + e ioo P = A f + A) where A = log(e ei °°p~ eii » k + 1). 

The free energy of an edge or a free face depends only 
on the effective thermodynamic potential of the polymers 
and the interaction between the stickers, which is similar 
in both the linked faces and free faces. We denote the 
free energy of an edge as OMfeilj), where f e {/j) is the 
free energy density (free energy per polymer site) of the 
face surface, or the free energy of one of the <tm sites 
on the face. Since in our model a free face differs from 
an edge only in the number of sites available (i.e., only 
ctm/2 polymers can be accommodated on a free face) and 
in the effective thermodynamic potential, its free energy 
is (cr M /2)/ e (> + A). 

One can readily derive the free energy, / e , of a system 
with only excluded volume interactions: f c = — l og(l + 
e~ M ). We considered such a system in Subsection IIII CI 
where we ignored the configurational free energy of the 
polymers. 

The interactions between the stickers are accounted 
for within a mean field approximation. The total inter- 
action energy for an edge with om polymers is written 
as —aaj\j/2. The number of stickers that occupy a given 
edge is a. We therefore write the interaction free energy 
due to two-body interactions among the stickers using a 
virial-type term: — (a/2aM)v 2 - The partition function is 



Z = J2 



cm , 

exp -^a 
a ) \ 2a m 



^ / d</> e-* ' 2a (e*^ + 1) (11) 
2ita 1 



We consider this expression for Z as the parti- 
tion function of an effective Hamiltonian, H(<f>) = 
au [^; < f >2 ~ log^e 1 ^ - ' 1 + 1)], that has a minimum at </>o, 
given by: 



a 



e ^-0o + 1 



(12) 



For aM 1 the distribution of <fi tends to a sharp peak 
around </>q, and we can use the mean field approximation 
(/(^))</> ~ f(4>o) to estimate the average fraction <p = 
(a)/aM of sites occupied in an edge and its variance: 



si 



df e = 1 / dH{4>) \ 

dfx a M \ dfi / ^ a 

d 2 f c ( 1 



dfi 2 



<p(l - ip) 



(13) 
(14) 
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FIG. 2: The variance of the number of polymers in a bond, 
as a function of the average bond occupation (determined by 
the chemical potential fi), for ctm = 20. The values of the 
attractive interaction, a, are in units of k^T. 



This expression applies whether a is positive (attractive 
interactions), negative (repulsive interactions), or zero. 
Eq. E| suggests a divergence of the variance for a > 4. 
In a macroscopic system (ctm ^ oo), this signifies an 
instability related to macroscopic phase separation. We 
do not expect such behavior for finite am (which is more 
realistic) , but we do observe a pronounced increase in the 
variance, as the interaction strength a increases (see Fig. 

The increase in the variance is due to a bimodal dis- 
tribution of the polymer density in the bonds, as demon- 
strated in Fig. |21 For large values of a, the distribution 
p(a) resembles the bimodal toy distribution in Eq. 0] 

In such a system, when the interactions are strong 
enough, one may observe two kinds of edges connecting 
the droplets: dilute edges containing one or a few poly- 
mers, and dense edges, containing about aM polymers. 
Both kinds induce the same effective attraction between 
the droplets (see below). Under degradation, most of the 
dilute edges will break down, and the connectivity of the 
network may decrease dramatically, as depicted in Fig. 
|TJl. If the system was in the gel phase, and the connec- 
tivity is reduced below the percolation threshold, the gel 
will transform into a fluid phase. 

We conclude by obtaining an analytical expression for 
the mean-field free energy of an edge, which is a good 
approximation for the case ctm S> 1. Expansion of the 
effective Hamiltonian around its minimum yields H(<f>) ss 
H(M + §£T"(0o)(0 - M 2 , where H"(<f>) = d 2 H/dtf. 
Plugging this approximation into the integral of Eq. 1111 
and taking the logarithm, we obtain the grand-canonical 
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FIG. 3: The distribution P(a) of the number of polymers in 
a bond, for om = 20 and <tm = 200. The effective chemical 
potential fj, is such that (a) = om/2. The values of a axe in 
units of ksT. 



potential for a single edge: 

/ e (A*) = -l0g(Z)w — 

= ^ 2 + log(l - 



#(0o) - r log f—H' 



Each edge has a free energy, <JMfe{n), and every free 
face has a free energy, (cr m /2) fc{f^ + A). The partition 
function of the system containing both polymers (which 
have both interactions and configurational entropy) and 
droplets (which have translational entropy) is given by 
the sum over all droplet configurations: 

Z = CX P [ N dVd + N K a M f e (jj) + N F ^f c (p + A 



Iog(l - a<p(l - <p)) , (15) 



2a m 



where tp is determined by the value of /i, and the last 
equality is obtained using Eg uat ions 1121 and 1131 



B. Droplet free energy 

Each vertex of the network in our model can oc- 
cupy one oil droplet. We define Sj = 0, 1 as the oc- 
cupation number associated with vertex i in the lat- 
tice. The number of edges in a given configuration of 
the system is = SjSj, where the sum is over 

pairs of neighboring sites. The number of free faces is 
N-p — Yltij) Si(l — Sj) + Sj(l — Si). The number of droplets 
is N d = £\ Si. 



^exp 



N d [fl d + ~<T M fe(V + A) ) + NEJiii, A) 



(16) 



where /id is the chemical potential of the droplets, 
J(/i, A) = f e (n) — /e(A* + A) is the effective interaction 
between the droplets, mediated by the polymers, and £ 
is the number of nearest neighbors at each site. This is 
equivalent to the partition function of a lattice gas model 



eq 
7] 



We use a Landau type, coarse grained approximation 
for the free energy, as a function of the dimensionless con- 
centration (volume fraction) of droplets, c = (s), which 
is the average occupation per vertex: 

/ = clogc+ (1 - c) log(l - c) 

+ ^C<Wc(M + A)c + ^(a M J(v, A)c 2 . (17) 

The first two terms in Eq. 1 171 give the lattice gas entropy; 
the third term is obtained using N d /V = c, where V is 
the volume of the system; the fourth term is obtained 
from Eq. I16l using the mean-field approximation N-e/V ~ 
(c 2 /2. The dependence of / on the potential /id is linear 
in the droplet concentration, and this term has no effect 
on the phase behavior of the system. 

In the coexistence region, there is an equilibrium of 
droplet chemical potential, [i d = df/dc, and of the os- 
motic pressure, 7r<j = / — fi d c, between the coexisting 
phases. In equilibrium, the polymer chemical potential 
must be the same in both phases; this is automatically 
ensured by working at constant /i. By solving these 
equalities we obtain the binodal curve in c x /z (droplet 
concentration-polymer thermodynamic potential) plane. 

Since in experiment it is easier to control the polymer 
concentration tp than the polymer chemical potential, it 
is useful to predict the phase diagram in the cx tp plane. 
We perform this transformation using the identity 



A) c +J'( M ,A)c 2 ] 



(18) 



where the prime indicates a derivative with respect to 
fi. This fixes the polymer chemical potential, /i, as a 
function of the number of polymers, p. In the phase- 
coexistence region, above the binodal curve, the condi- 
tions of constant droplet chemical potential, /id, and os- 
motic pressure, 7Td, for fixed polymer chemical potential, 
/i, yield values of the droplet density, c, for the coexist- 
ing dilute and dense phases respecitvely. From Eq. 1181 
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FIG. 4: The phase coexistence curves of systems with differ- 
ent values for the interaction between the stickers, a (given 
in units of ksT). The curves are obtained from the condition 
of equilibrium of the droplet chemical potentials, fid, and os- 
motic pressure, 7Td, in the coexisting phases, using the Lan- 
dau free energy is given by Eq. 1171 The a;-axis represents the 
droplet volume fraction, c, and the y-axis represents the num- 
ber of polymers per available site, (p. The parameters used 
are A = 0.3fc B r, £ = 4 and a M = 20. 



this implies that each phase has a different value of the 
polymer concentration, tp. 

The critical point in the c X fi (droplet concentration- 
polymer concentration) plane occurs at c = 1/2, with fi 
given by 



J(/i,A) 



(19) 



The binodal curve of the system is presented in Fig. 01 
for different values of the edge interaction a. The max- 
imal number of stickers that can reside in one droplet 
is Qijm- Thus, the number of stickers per unit volume is 
cQctm , and the maximal number of polymers per unit vol- 
ume is c(ctm/2 (assuming a prohibitively high energetic 
cost of a hydrophobic sticker in the water). The number 
of polymers per available site in the system is then 



2(p 
(cr M c 



(20) 



A phase co-existence region is present in a system with- 
out explicit sticker interactions (a = 0) 0). However, the 
interactions can markedly affect the shape of the coex- 
istence curve and of the concentrations in the two co- 
existing phases. As the interaction, a, increases, the ef- 
fective interaction between the droplets, induced by the 
edges, increases and becomes more sensitive to the poly- 
mer concentration. As a result, the onset of phase sepa- 
ration is sharper, and characterized by a phase with very 
dilute droplets and low polymer concentration that co- 
exists with a phase with high droplet density and high 



polymer concentration. In fact, in Fig. 0] the curve for 
a = 10 is almost parallel to the tie-lines. 

We note that dilute edges (with one or a few polymers) 
can also exist in the dense droplet phase, since they in- 
duce the same effective interaction, J(fi, A), between the 
droplets as the dense edges, and their polymers have the 
same thermodynamic potential, /i, as the polymers in 
the dense edges. This scenario may change if we consider 
interactions between the edges. 



V. DISCUSSION 

In this paper, we have predicted the immunity of the 
gel phase of a polymer-microemulsion system to random 
degradation of polymers. The gel in this system consists 
of oil droplets connected by tclcchclic polymers. The set 
of polymers that link a given pair of droplets is defined 
as an edge (i.e., bond) of the gel network. We consider 
the distribution, p(cr), of the number of polymers in one 
such edge, and argue that immunity of the gel is inversely 
related to the variance, S%, of this distribution. 

Repulsive interactions between the stickers (the 
tclcchelic ends of the polymers) reduce the variance S%, 
thus contributing to the gel immunity. Attractive inter- 
action, which increase S%, reduce the immunity of the 
gel. 

Attractive interactions may give rise to a bimodal dis- 
tribution of the edge occupation, which can give rise to 
the co-existence of dense and dilute edges in the same 
macroscopic sample. A dilute edge contains a small num- 
ber of polymers, and is likely to vanish under a random 
degradation of the polymers. When a finite fraction (say, 
1/2) of the edges are dilute, the fraction of edges surviv- 
ing a degradation will be lower, and the system is more 
likely to transform from the gel to the fluid phase. 

This reasoning is particularly applicable to the case 
of a dense phase that coexists with a dilute phase (as- 
suming the dilute phase is not a gel). Since the polymer 
density within the dense droplet phase is high, we ex- 
pect - if the attractive interactions are large enough - to 
find a bimodal distribution of the edge occupation within 
the dense droplet phase. In that case, we are interested 
in the immunity of the dense phase, and therefore must 
take into account the distribution, p(a), within the dense 
droplet phase. 

To test our predictions for the immunity of a system 
with attractive interactions among the polymers, we sug- 
gest comparing the elastic stability (gel-like nature) with 
respect to polymer degradation of two systems (both in 
the dense phase): one with no sticker interactions (a = 0) 
and one with strong interactions (say, a = 10). Accord- 
ing to our predictions, the system with no interactions 
will be characterized by edges each of which has approx- 
imately the same number of polymers, while the system 
with strong interactions will show a bimodal distribution 
for the edge occupancy. We predict that after polymer 
degradation the fraction of surviving edges will be larger 
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in the non-interacting system (see Fig. and that this 
system is more likely to remain in the gel phase. 

In order to investigate the immunity of such gels exper- 
imentally it is necessary to devise a method to degrade 
the polymers at random and in a homogeneous manner 
throughout the system. It may be difficult to do this us- 
ing molecules in a solvent, since a solvent added to the 
system will initially have a much higher concentration at 
the surface of the gel (penetrating the bulk by slow diffu- 
sion); thus the polymers near the surface are much more 
likely to be degraded. 

One way to achieve random degradation may be to 
use polymers that have been synthesized with a photo- 
sensitive molecule in the middle of each chain. Opti- 



cal methods might be used to cause these molecules to 
change their conformation and dcpolymcrizc. The in- 
tensity and duration of the optical beam can be used to 
control the fraction of polymers broken. 
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